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§1.  Introduction 

T3m  pocking  of  spheres  has  been  a  subject  of  interest  far  many  ■«—**«*«  (Sloane 
(1084)).  Many  investigators  hare  sought  the  densest  way  to  arrange  identical  spheres  in 
space.  An  asymptotic  formula  that  provides  an  upper  bound  for  the  packing  density,  94, 
of  dimensional  spheres  in  a  dimensional  space  aa  the  space  volume  approaches  infinity 
is  doe  to  HJS.  Daniels  and  reported  in  Rogers  (1988, 1964).  It  is 


In  addition  several  investigators  have  looked  into  the  random  packing  of  spheres.  Ran¬ 
dom  sequential  packing  of  spheres  hae  been  appSed  by  Bernal  (1988)  to  study  the  structnxe 
of  liquids  and  since  then  it  has  been  rilsmsssri  by  other  authors,  for  example,  HIguti  (1960), 
Solomon  (1987),  Thaemnra  (1979)  and  others.  Because  of  analytical  difficulties  when 
a  >  2,  one  dimensional  random  sequential  packing  has  received  attention.  It  is  an  inter¬ 
esting  subject  in  probability  and  has  been  dbcureed  by  several  authors,  for  example,  Rfinyi 
(1988),  Dvoretsky  and  Robbins  (1964),  Solomon  (1967),  and  Rah  (1980).  The  one  dimen¬ 
sional  model  can  be  extended  to  random  sequential  of  cubes.  Consider  random 

sequential  paddng  of  hypercubes  ci  ridabngth  1  brio  a  larger  hypercube  of  akfabngth  s 
with  rigid  boundaries.  R  seems  to  be  natural  to  expect  that  the  B"»«ti"g  density  As  aa 
s  tends  to  infinity  in  a  space  of  dlmenriom  i  squab  0*  where  0  is  the  limiting  packing 
density  for  i  »  L 

This  conjecture  was  poeed  and  dboieiad  by  Palasti  (1980).  as  a  follow  up  to  Rinyi’a 
pioneering  artide.and  his  results  for  one  dimensional  sequential  random  paddng.  Several 
emulations  and  analyses  by  Blabddl  and  Solomon  (1970, 1982)  and  Akeda  and  Hoti  (1975, 
1978)  refute  her  conjecture.  Asymptotic  behavior  of  the  random  paddng  density  aa  the 
dimension  tends  to  infinity  b  of  intend,  but  thee  b  no  analytical  result  for  dimensions 
higher  than  one.  Computer  snmilatoos  are  mS  expensive  as  increases. 

Here  we  introduce  two  sample  models:  shnpb  cube  random  paddng  and  random  paddng 
ky  Hamming  distance,  to  get  some  feel  for  random  parking  density  hi  higher  dimension. 
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We  improve  the  previous  «mniation  results  for  cable  random  packing  (Bob  and  Ueda 
(1983))  by  increasing  tbs  number  of  runs.  Also,  we  give  aimnlatfon  results  for  random 
packing  by  distance  and  dismiss  the  behavior  of  packing  density  when  dimen¬ 

sionality  is  increased.  Ibr  the  earn  of  Hamming  distances  of  2  or  3,  d“*  fits  the  rimnlatfoa 
results  of  packing  density  where  d  represents  the  dimension  and  is  an  empirical 
constant  astiaatad  by  tha  laast  squares  method, a aa  Table  3. 

The  variance  of  packing  density  is  larger  when  fe  is  even  and  «"**n**‘  when  k  is 
odd,  see  Table  2,  where  A  represents  Hamming  distance . 

§2.  Simple  enble  random  sequential  pecking 

BUWdI  Soto*  (MB)  |m  »  npoimatil  fonacla 

tf'-fimid-l Wp-fil)  far 

dam  3, 4,  which  may  direct  oer  attention  to  packing  density  bi  higher  dimensions.  Bet  for 
the  more  general  stedka  dlscnssed  by  the  above  anthem.  It  is  dUBcuit  to  get  vainas  for 
mw  ttmi  fonr  dimensions.  For  a  simple  random  packing  ya^h  «««i  Uiada  (1083) 

gave  parking  densities  ep  to  11  dlmanrieaa.  In  this  model  consider  a  cabs  in  d  dimensions 
with  eideiength  4  end  a  cubic  lattice  with  unit  aiddength.  Cubes  of  ridstagth  2  are  pvt 
sequentially  at  random  into  the  cobe  of  aiddength  4  so  that  each  vertex  coincides  with 
one  of  the  lattice  points.  We  ccntfame  until  no  place  can  be  found  in  the  large  cube  to 
place  the  smaller  one.  Consider  the  packing  density  fg  of  dimension  d.  The  experimental 
formula  by  Blaisddl  and  Solomon  (1982)  fits  the  results  Sard  •  3, 4,5  aa  we  can  eaaOyaee 
from  the  results  by  Boh  and  Ueda  (1983).  SVom  computer  nmnlationa  up  to  dimension 
11,  Tsm/ts  seams  to  approach  one  (Boh  and  Ueda  (1983)),  which  may  suggest  04+ilfk 
of  the  general  modal  also  approaches  one  as  d  tends  to  infinity,  see  Thbla  L 

$3.  Random  sequential  coding  by  Hemming  distance 

One  of  the  poarible  applications  of  spherical  random  packing  la  in  recognition  theory, 
(Dolby  and  Solomon  (1978)).  To  mtutmi—  the  recognition  error  rate  for  a  given  rejection 


cate,  the  following  decision  rale  is  reasonable,  if  the  objects  axe  all  equally  Hbdy  to  occnr. 
If  an  arbitrary  point  in  the  observation  space  is  within  r  distance  units  of  the  expected 
value  point  of  one  (and  only  one)  of  the  given  objects,  identify  the  arhitraxy  point  with  that 
object.  Otherwise  reject  the  observation.  Here  distance  is  the  nsnal  Euclidean  distance. 
Dolby  and  Solomon  (1075)  compared  spherical  random  packing  density  with  the  ratio 
between  existing  monosyllabic  words  and  possible  monosyllabic  forms  in  Here  we 

simplify  the  situation  and  introduce  a  random  sequential  coding  model. 

Consider  a  set  of  2*  points  whose  coordinates  are  1  or  0  in  a  Euclidean  space  of 
dimension  d.  Euclidean  distance  Is  defined  between  two  points  of  the  2*  points.  The 
square  of  the  Euclidean  distance  in  this  case  is  called  the  ffawMwfag  distance  in  coding 
theory.  Problems  in  random  coding  axe  discussed  by  several  authors,  (Shannon  (1948), 
Abramson  (1983)  ).  Our  model  may  be  called  random  sequential  coding.  Consider  a 
random  sequential  packing  into  the  2*  points.  At  first  we  chooee  one  point  (d  coordinates) 
at  random  and  we  record  it.  Chooee  another  and  record  it  if  its  distance  is 

^  h,(h  <  d)t  otherwise  discard  it.  How,  choose  the  next  point  at  random  and  record  it 
if  the  Hamming  distance  from  each  of  the  2  points  is  not  less  than  h,  otherwise  discard  it 
and  chooee  another  point  at  random.  We  continue  this  procedure  until  there  is  no  possible 
point  to  record  among  the  2*  points  and  we  now  have  the  number  of  recorded  points 
X(d,Jk). 

Define  a  packing  density  by  X(d,Jfc)/24  and  label  X(X(d,A)/ 2*)  *  64*.  Computer 
simulations  suggest  64+  approaches  one  as  d  tends  to  infinity,  as  in  the  case  of  the 

simplest  cubic  random  padring.  Another  interesting  fact  is  the  variation  of  the  variance  of 
X(n,k).  For  k  even  the  variance  is  large  and  for  odd  k  the  variance  is  small;  for  fixed  d. 
sat  Table  2. 

Our  model  will  help  to  the  coding  system  arising  in  nature  (as  well  as  finguie* 

tie  problems).  Ear  example,  in  the  case  of  the  ammo  add  code,  64  words  are  theoretically 
possible  in  the  triplet  coding  system  by  four  species  of  nucleotides  La.  4*.  The  actual 
number  of  amino  adds  plus  chain  terminator  in  the  code  is  only  21.  The  ratio  of  the 
number  of  elements  actually  used  to  the  number  possible  is  21/64  •  0X2812$.  ft  should 
be  noted  that  this  value  is  quite  close  to  0X283,  the  paddsg  value  for  d  *  6  and  km  2. 
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Note  that  2*  »  64.  A  random  sequential  packing  model  by  Dolby  and  Solomon  (1975) 
waa  applied  to  this  problem  by  Bah  and  Hasegawa  (1980).  But  onr  present  model  may  be 
more  amenable  far  this  situation. 


$4.  Packing  density  sad  experimental  formula 
Rom  the  inequality  far  Bamming  bounds,  we  hare 


far  k  m  2s  + 1,  see  Lin  (1988).  Hence,  if  the  limit  of  exists,  it  shoold  be  one. 

The  Hamming  upper  bound  can  be  strengthened  aa  fallows.  Let  M{dtk)  be  the  w«winm 
possible  number  at  code  words  in  a  binary  code  with  words  at  length  i  and  «*» tiwi» 
distance  k,  if  k/d  <  then 


M^k)  Z 


d  l*JT(p) 
22?  ©(£-') 


when  p  m  i/2 d,i  as  (l/2p(l  -  Vl~4pj,JT(p)  *  p/>/ l-4p  (Wyner  (1984),  Ash  (1989)). 
Consider  a  regular  packing  far  k  m  2.  Consider  a  set  of  points  on  the  hyperpianee  st  + 


*»  +  •••**  »  m,  where  m  is  odd.  Then  the  mutual  distance  between  each  pair  of  points 
on  the  set  is  not  lass  than  2.  Thbe  every  point  on  the  hyperplaae  far  odd  m  and  do  not 
take  any  point  at  the  byperplanae  far  even  m.  This  makes  a  regular  packing  of  density 
OS,  which  does  not  depend  on  dhnenihni.  The  farm  d“*  fits  oar  random  paridng  density 
reasonably  well  for  h  ^  2, 3,  as  we  shall  see  in  Table  3,^<r«  «  i*  estimated  fro*  the 


logarithmic  plot  of  the  simulation  result  for  10^<f<17  by  tha  laast  squares 
msthod.  For  the  case  Jk  >  4,  aT*  does  not  fit  our  simnlation  results. 


We  would  like  to  thank  Joe  Marhoul  and  Akin  Yoehida  for  their  hdp  in  programming 
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Sequential  random  packing.  Cubic  packing,  Coding  theory.  ^ 
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J  Simple  cubic  random  packing  and  ranaom  packing  by  jamming  distance  are  ana- 
|  lyzed.  Consider  the  packing  density  y#  of  dimension  by  cubic  random  packing. 

I  From  simulations  up  to  dimension . 11,  /y^f seems  to  approach  one.  Simulation 

results  for  random  packing  by  Hamming  distance  are  obtained.  For  the  case  of 
Hamming  distances  of  2  or  3,  fits  the  simulation  results  of  packing  density 
where  a  is  an  empirical  constant. 
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